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MAGNETIC  STUDIES  OF  THE  ANTIFERROMAGNET  RbMnF3 


ABSTRACT 


Magnetic  properties  of  the  antiferromagnet  RbMnFg  have  been  studied.  Below  TN' 
the  magnetic  ions  order  into  a  two-sublattice  system  with  the  spins  antiparallel . 
RbMnFg  exhibits  high  exchange  and  low  anisotropy;  the  form  of  the  anisotropy  sur¬ 
face  is  cubic.  Consequently,  for  applied  DC  magnetic  fields  less  than  about 
3000  oe,  the  static  equilibrium  position  of  the  sublattice  magnetization  is,  in 
general,  multivalued. 


Measurements  of  the  DC  susceptibility  y  have  been  made  on  powder  and  single 
crystal  specimens  of  RbMnF^  for  the  range  of  applied  field  0  to  12  koe  and  over  the 
temperature  range  4.2  to  300°K.  The  observed  value  of  T^  was  about  ten  degrees 
higher  than  the  previously  published  value.  When  plotted  as  a  function  of  applied 
field,  Xjiiij  shows  no  abrupt  discontinuity  analogous  to  the  spin  flopping  exhibited 
by  uniaxial  antiferromagnets. 


A  simple  model,  in  which  H^  and  M  are  restricted  to  the  (110)  plane,  has  en¬ 
abled  solutions  of  the  static  equilibrium  problem  to  be  obtained.  X-band  resonance 
experiments  are  reported,  and  a  resonance  theory  is  presented  which  incorporates 
the  equilibrium  solutions.  The  predicted  antiferromagnetic  resonance  spectrum 
shows  reasonable  agreement  with  the  experimental  data.  The  possibility  of  parallel 
pumping  spin  waves  in  RbMnFg  is  considered,  and  an  attempt  to  measure  the  spin 
wave  linewidth  is  described. 
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MAGNETIC  STUDIES  OF  THE  ANTIFERROMAGNET  RbMnFg 

CHAPTER  I 
INTRODUCTION 

This  report  describes  certain  aspects  of  the  static  and  dynamic  properties  of  RbMnFg  —  a 

perovskite  material  which  has  a  simple  cubic  structure  and  is  antiferromagnetic  at  temperatures 

2+ 

less  than  82.6°K.  Below  the  Neel  temperature,  the  Mn  ions  order  into  a  two-sublattice  con¬ 
figuration  with  the  magnetic  moments  antiparallel.  According  to  published  data, 1,2  RbMnF- 
exhibits  low  anisotropy  and  high  exchange  fields  [356  and  7.08  X  10  amp/meter  (4.47  and 
8.9  X  105oe),  respectively!.  Thus,  antiferromagnetic  resonance  (AFMR)  experiments  can  be 
performed  without  the  need  for  extremely  high  DC  magnetic  fields  and  at  convenient  microwave 
frequencies  (e.g.,  X-band).  The  form  of  the  anisotropy  appears  to  be  entirely  cubic,  with  no 
uniaxial  component.  The  anisotropy  constant  is  negative,  so  that  the  four  cube  diagonals  are 
equivalent  easy  directions  of  magnetization. 

Because  RbMnFg  is  chemically  quite  stable,  it  is  a  very  convenient  material  for  general 
antiferromagnetic  studies.  It  is  possible  to  prepare  relatively  large  single  crystals  which  may 
be  readily  cut  by  conventional  techniques.  In  the  single  crystal  form,  the  material  shows  little 
or  no  tendency  to  absorb  moisture. 

The  equilibrium  position  of  the  sublattice  magnetizations  under  the  action  of  an  applied  DC 
magnetic  field,  determined  by  minimizing  the  total  free  energy  of  the  system,  is  not,  in  general, 
unique  for  an  antiferromagnet  that  possesses  cubic  anisotropy.  However,  for  sufficiently  large 
values,  the  applied  field  and  the  magnetization  will  be  almost  orthogonal,  regardless  of  the 
direction  of  the  applied  field  with  respect  to  the  cube  edges.  In  this  case,  the  equilibrium  posi¬ 
tion  of  the  sublattice  magnetizations  is  well  defined.  For  RbMnFg,  experimental  evidence  shows 

that  the  perpendicular  configuration  is  quite  an  accurate  model  for  fields  in  excess  of  8  koe. 
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Freiser,  et  al.,  have  performed  a  comprehensive  investigation  of  AFMR  for  the  high  field  case 
which  has  enabled  them  to  investigate  the  anisotropy  surface. 

To  the  author's  knowledge,  no  work  has  been  published  concerning  the  nature  of  the  AFMR 

„  3 

spectrum  in  RbMnFg  in  the  low  field  region.  For  weak  fields,  Neel  has  postulated  the  existence 
of  domains  in  cubic  antiferromagnets  which,  since  the  dipolar  fields  are  small,  are  probably  ill 
defined.  Because  of  the  likelihood  of  their  existence,  the  resonance  spectrum  can  be  expected 
to  be  much  more  complicated. 

Measurement  of  the  DC  magnetic  susceptibility  as  a  function  of  applied  field  for  various 
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crystal  axes  gives  some  insight  into  the  equilibrium  position  of  the  magnetization.  McGuire 
has  determined  Xj^qO]  at  anc*  *  k°e  by  measuring  the  force  on  a  single  crystal  of  RbMnFg 


placed  in  a  nonuniform  magnetic  field.  He  found  X j 1 00]  to  °-8X{  and  concluded  that  the  spins 
were  distributed  equally  along  the  four  cube  diagonals.  This  simple  picture  cannot  be  correct, 
since  the  torque  due  to  the  applied  field  must,  at  least,  rotate  the  magnetization  vectors  toward 
the  [110]  axis. 

A  study  of  the  DC  magnetic  susceptibility  was  made  using  a  vibrating  sample  magnetometer. 
The  results  are  presented  in  Ch.  II.  Measurements  at  4.2°K  were  made  on  a  powder  sample  as 
well  as  on  a  single  crystal.  The  results  show  that  Xj^qqj  and  Xjmj  are  both  increasing  functions 
of  but  approach  an  asymptotic  value  at  about  8koe.  No  effect  was  observed  which  corre¬ 

sponded  to  the  abrupt  spin  flop  transition  of  uniaxial  antiferromagnets.^ 

Susceptibility  measurements  were  also  recorded  as  a  function  of  temperature,  from  4.2°  to 
300°K,  enabling  the  Neel  temperature  and  the  exchange  constants  and  to  be  determined. 

Despite  the  complex  nature  of  the  problem  concerning  the  static  equilibrium  position  of  the 
magnetization,  the  situation  is  considerably  simplified  if  is  restricted  to  lie  in  one  of  the 

crystal  symmetry  planes.  A  discussion  of  the  static  equilibrium  is  presented  in  Ch.  III.  The 
analysis,  based  on  the  molecular  field  approximation,  assumes  that  the  magnetization  is  uniform, 
i.e.,  the  discussion  is  limited  to  a  single  domain,  and  that  H^^,  lies  in  a  (110)  plane.  As  ex¬ 
pected,  even  for  the  restricted  problem,  the  solution  is  multivalued. 

A  discussion  of  AFMR  for  a  cubic  material  is  presented  in  Ch.  IV.  Analysis  of  the  linear¬ 
ized  equations  of  motion  yields  two  independent  solutions  for  the  resonant  frequency  of  the  system. 
The  particular  solutions  obtained  when  H^^,  is  chosen  to  lie  along  a  [111]  axis  are  derived  from 
an  equation  which  is  equivalent  to  the  result  obtained  by  Keffer  and  Kittel f* 

X-band  resonance  data  were  obtained  at  five  different  frequencies.  The  resonances  were 
studied  in  detail  at  4.2°K  as  a  function  of  angle  and  Hj^.  Plots. of  Hj-^  required  to  excite  reso¬ 
nance  vs  angle  of  applied  field  with  respect  to  the  cube  edges  indicate  at  least  two  resonances  for 
each  frequency.  For  at  least  two  frequencies,  many  resonances  are  apparent  which  are  clearly 
inexplicable  on  the  basis  of  the  simple  theory.  However,  at  least  one  resonant  mode  has  been 
correlated  as  a  function  of  frequency.  Some  of  the  characteristics  of  this  mode  can  be  predicted 
and  are  in  agreement  with  the  theory  presented. 

Spin  wave  instability  and  premature  saturation  of  AFMR  have  been  discussed  by  Heeger  and 
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Pincus.  The  mechanism  of  spin  wave  growth  and  instability  was  first  analyzed  for  the  ferro- 

g 

magnetic  case  by  Suhl.  He  showed  that  such  instabilities  could  be  explained  by  a  coupling  of 
energy  from  the  uniform  precessional  mode  (spin  wave  number  k  =  0)  to  spin  waves  (k  7^0)  which 
are  degenerate  in  energy  with  the  uniform  mode.  The  spin  waves  in  turn  couple  energy  to  dissi¬ 
pative  processes  such  as  lattice  vibrations.  When  the  threshold  driving  field  for  the  excitation 
of  spin  waves  has  been  exceeded,  the  amplitude  of  the  uniform  precessional  mode  remains  con¬ 
stant.  Any  further  increase  in  incident  power  is  coupled  to  the  spin  wave  spectrum.  The  analy¬ 
sis  by  Heeger  and  Pincus  for  the  antiferromagnetic  case  is  similar.  By  considering  the  equations 
of  motion  for  a  two-sublattice  system,  they  have  calculated  the  critical  field  for  the  onset  of 
instability  and  premature  saturation  of  the  AFMR. 


The  spin  wave  linewidth  in  KMnF^,  a  material  similar  in  structure  to 


RbMnF. 


has  been 


measured  by  Heeger.  He  observed  the  apparent  saturation  of  the  susceptibility  of  the  uniform 
mode  at  resonance  and  demonstrated  that  the  critical  field  was  anomalously  low.  The  experi¬ 
mental  data  indicated  that  AH^  was  about  5X10  oe,  which  is  nearly  four  orders  of  magnitude 
less  than  the  linewidth  of  the  uniform  precessional  mode;  AH^  was  inferred  from  a  measurement 
of  the  response  of  the  resonant  sample  to  amplitude  modulated  incident  power. 
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Heeger's  measurements  indicate  that  the  mechanisms  responsible  for  the  line  broadening  in 
spin  wave  excitation  must  be  very  different  from  those  associated  with  the  uniform  mode  for  the 
material  KMnF^.  The  question  arises  whether  RbMnF^  exhibits  similar  behavior  and  whether 
the  spin  wave  linewidth  is  the  same  order  of  magnitude.  Chapter  V  describes  a  method  and  an 
experimental  attempt  to  obtain  AH^  in  RbMnF^  by  measuring  the  threshold  magnetic  field  for 
longitudinally  pumped  spin  waves.  The  attractive  feature  of  the  so-called  parallel  pumping 
technique  is  that  the  uniform  precessional  mode  is  suppressed.  This  permits  a  more  accurate 
measurement  of  spin  wave  linewidth,  since  the  losses  associated  with  the  uniform  precessional 
mode  are  absent.  The  result  of  the  experiment  was  negative,  which,  however,  was  not  inter¬ 
preted  as  having  set  a  lower  limit  to  the  value  of  AH^.  Analysis  of  the  instability  threshold  makes 
the  basic  assumption  that  the  pumping  field  is  parallel  to  the  magnetization  and  that  the  spin  wave 
frequencies  excited  are  very  low.  The  negative  result  of  this  experiment,  together  with  the  re¬ 
sults  of  the  DC  susceptibility  and  AFMR  experiments,  was  taken  to  imply  that  these  conditions 
were  violated. 
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CHAPTER  II 

DC  SUSCEPTIBILITY  MEASUREMENTS 


A.  GENERAL  DISCUSSION 


The  sublattice  magnetizations  of  RbMnF^  are  quite  well  described  by  the  Brillouin  function 

for  J  =  5/2.  The  saturation  magnetization  M~  has  been  calculated  from  the  formula,  assuming 

^  o  *  5  , 

a  g-factor  of  two  and  the  lattice  constant  equal  to  4.24  A,  and  was  found  to  be  3.04  X  10  amp/ 

meter  (4nM^  equal  to  3820  gauss). 

According  to  the  molecular  field  theory  as  applied  to  antiferromagnets,  the  fields  acting  on 
the  respective  sublattices  are 


where 


H . 


H, 


=  H 


DC 


=  H 


DC 


+  W,  ,M 
1 1 

+  W22M 


1 

2 


+  W12M 
+  W0.M 

c  1 


2 

1 


W  =  W 
W11  22 


and 


W  =  W 
w  1 2  W21 


The  sum  of  the  exchange  constants  and  W^,  defined  by  the  above  equation,  may  be 

evaluated  from  a  plot  of  l/x  vs  temperature  within  the  paramagnetic  region.  The  plot  is  a 
straight  line  of  the  form 

7  =  T<25>-i<w„*w12> 


where 


2 

eff'  ’  o 


n(P^J  p 


3k 


(1) 


2+ 

k  is  Boltzmann’s  constant,  n  is  the  number  of  Mn  ions  per  unit  volume,  and  P  ff  is  the  mag- 

24-  en 

netic  moment  of  the  Mn  ion. 

The  straight  line  plot  intercepts  the  negative  temperature  axis  at  a  point  given  by 


0  =  C(WU  +  W12) 


(2) 


The  value  of  the  slope  and  intercept  enable  (W^  +  W12)  to  be  found;  and  may  be  com¬ 

pletely  determined  from  a  measurement  of  the  large  field  susceptibility,  which  is  equal  to  l/W.7, 

10 

below  T The  magnitude  of  the  slope  of  the  line  enables  P  whose  value  is  well  known,  to 
be  verified.  This  provides  one  check  of  the  accuracy  of  the  measurements. 

In  addition  to  providing  values  of  T^,  and  W measurements  of  the  DC  susceptibility 

can  give  some  insight  into  the  equilibrium  position  of  the  magnetization. 


*  The  quantum  defect  has  not  been  taken  into  account. 
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A  series  of  measurements  of  the  magnetic  susceptibility  x  was  made  on  powder  and  single 
crystal  specimens  of  RbMnF^.  The  single  crystal  was  cut  so  that  it  could  be  rotated  about  an 
axis  perpendicular  to  the  (110)  plane.  Measurements  of  x  with  Hj-^  applied  parallel  to  the  [111] 
and  [100]  axes  enabled  two  determinations  of  ^  and  W^.  In  the  case  of  a  powder,  the  density, 
which  enters  into  the  calculation  for  the  exchange  constants,  is  not  well  defined.  Hence,  the 
powder  measurements  were  not  used  to  determine  and  W^. 

B.  EXPERIMENTAL  ARRANGEMENT 

Magnetic  susceptibility  measurements  were  made  using  a  vibrating  sample  magnetometer 

11  12 

previously  described  by  Hunt  and  more  recently  modified  by  Santoro.  A  block  diagram  of 

13 

the  instrument,  which  is  modeled  after  Foner's  design,  is  shown  in  Fig.  1.  The  magnetometer 
measures  the  magnetic  moment  of  a  sample  placed  in  a  uniform  magnetic  field.  The  sample  is 
made  to  vibrate  in  a  direction  perpendicular  to  the  DC  field  by  a  voltage  which  is  induced  in  a 
pair  of  stationary  pickup  coils  that  surround  the  sample.  The  vibrational  motion  of  the  sample 


| 3-44-8867  I 


Fig.  1 .  Block  diagram  of  magnetometer. 


*  All  the  RbMnF3  samples  used  in  this  work  were  grown  and  donated  by  Dr.  M.  Kestigian  of  the  Sperry  Rand 
Research  Center,  Sudbury,  Mass. 
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Fig. 2.  x/p  vs  Hq£  for  RbMnF3 
powder  specimen. 
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originates  in  a  loud  speaker  cone.  The  sample  is  fastened  to  the  end  of  a  rod  which  is  in  turn 
rigidly  attached  to  the  cone.  A  second  pair  of  pickup  coils  surrounds  the  upper  part  of  the  sup¬ 
port  rod.  A  small  piece  of  ferrite  attached  to  the  rod  at  this  point  induces  a  reference  voltage 
in  the  second  pair.  By  the  use  of  amplitude  and  phase  adjustments,  the  signal  and  reference 
signals  are  balanced  against  each  other,  after  appropriate  amplification,  to  obtain  a  null  condition. 
The  magnetization  per  unit  mass  of  the  sample  is  then  equal  to  RK/W,  where  K  is  a  constant  re¬ 
lated  only  to  the  magnetometer,  W  is  the  weight  of  the  sample,  and  R  is  the  attenuator  reading. 
The  factor  K  is  found  experimentally  by  calibrating  the  instrument  with  a  sample  whose  magneti¬ 
zation  is  known.  The  actual  standard  used  was  a  nickel  sphere  (x/p  =  54.4emu/gm  at  300°K  and 
is  independent  of  applied  field  for  >  2  koe). 

The  constant  K  is  a  function  of  sample  position  and  the  appropriate  electromagnetic  con¬ 
version  factors.  As  a  function  of  position,  K(x,  y,  z)  has  a  saddle  point.  It  is  necessary  to  locate 
the  sample  as  near  as  possible  to  this  point;  otherwise,  a  slight  displacement  of  the  sample  causes 
a  change  in  the  value  of  K.  This  is  an  important  consideration  in  order  to  obtain  good  reproduci¬ 
bility  of  measurements  and  accuracy  of  calibration. 

For  the  purpose  of  making  susceptibility  measurements  as  a  function  of  temperature,  the 
sample  is  cooled  to  liquid  helium  temperature  and  then  allowed  to  return  to  room  temperature 
over  a  period  of  several  hours.  Since  there  is  little  room  between  the  pole  pieces  of  the  magnet, 
the  magnetometer  is  equipped  with  a  special  narrow-necked  dewar  of  somewhat  limited  capacity. 
The  sample  rod  extends  down  into  the  dewar  neck,  and  the  helium  is  drawn  up  around  the  sample 
through  a  glass  capillary  tube  that  encases  the  sample  rod  and  extends  down  into  the  body  of  the 
dewar  vessel.  The  helium  supply  can  maintain  the  sample  at  4.2°K  for  about  10  minutes. 

Temperature  is  measured  by  means  of  a  copper  constantan  thermocouple.  For  the  lowest 
temperatures,  the  liquid  helium  itself  is  used  as  the  reference  bath.  As  the  sample  warms  and 
the  helium  in  the  reservoir  becomes  depleted,  an  ice  bath  is  adopted  as  the  reference. 

C.  EXPERIMENTAL  RESULTS 

The  magnetic  susceptibility  per  unit  mass  x/p  (p  =  density)  of  a  powder  specimen  was  meas¬ 
ured  as  a  function  of  DC  magnetic  field  at  4.2°K  over  the  range  1  to  11  koe.  The  results  are 
plotted  in  Fig.  2.  It  can  be  seen  that  x/p  rises  rapidly  at  first  and  then  levels  off  at  about  3  koe 
Above  3  koe,  x/p  increases  almost  linearly,  but  at  a  much  slower  rate. 

The  susceptibility  at  10  koe  was  measured  as  a  function  of  temperature  over  the  range  4.2° 

to  300°K.  The  results,  shown  in  Fig.  3,  clearly  illustrate  the  Neel  temperature.  The  suscepti- 

14 

bility  in  the  antiferromagnetic  region  follows  the  Van  Vleck  theory.  This  indicates  that  although 
the  sublattice  magnetizations  are  antiparallel  within  each  domain,  the  static  equilibrium  position 
of  the  magnetization  vectors  must  be  randomly  distributed  throughout  the  powder.  The  single 
crystal  anisotropy  is  ineffective  in  determining  the  orientation  of  the  sublattices  with  respect  to 
^DC’  ^  ^  were  effective,  the  spins  would  assume  a  position  almost  perpendicular  to  Hj-^  for  a 
field  as  large  as  10  koe,  and  x/p  would  be  fairly  independent  of  temperature. 

The  considerable  discrepancy  between  the  measured  value  of  T^  of  94°K  and  the  previously 
published  value  of  82.6°K  was  at  first  attributed  to  the  fact  that  the  thermocouple  junction  was 
not  in  direct  contact  with  the  powder  sample.  This  source  of  error  was  later  discounted  because 
the  same  discrepancy  was  noted  in  the  single  crystal  experiments  in  which  it  was  possible  to 
locate  the  junction  much  closer  to  the  sample. 
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A  plot  of  (x/p)  1  vs  temperature,  shown  in  Fig.  4,  yielded  a  straight  line  in  the  paramagnetic 

region.  From  the  slope  of  the  line  and  Eq.  (1),  P  was  calculated  to  be  5.96  Bohr  magnetons, 

en  10 

which  is  in  excellent  agreement  with  the  value  of  5.92  quoted  by  Bates. 

The  powder  specimen,  which  weighed  0.82  gm,  was  prepared  by  grinding  down  a  quantity  of 
polycrystalline  material  and  was  sealed  in  a  phenolic  holder  to  prevent  water  absorption.  Be¬ 
cause  of  its  size,  not  all  the  sample  contributed  fully  to  the  magnetic  moment.  A  correction 
factor  was  applied  which  was  equivalent  to  using  an  effective  mass,  rather  than  the  actual  mass, 
in  the  susceptibility  calculations.  The  correction  factor  was  determined  from  a  comparison  with 
susceptibility  data  obtained  from  experiments  with  a  single  crystal  sample  which  was  much 
smaller  (0.48  gm). 

Figure  5  shows  the  variation  of  x/p  with  H for  a  single  crystal  oriented  so  that  H^^  was 
applied  along  the  [111]  axis.  The  figure  indicates  a  rapid  increase  in  x/p  as  H^^,  is  increased 
from  1  to  4  koe,  followed  by  a  gradual  asymptotic  leveling  off  to  a  value  of  8.55  X  10  ^  emu/gm. 


=  47r  X  Density  X  (x  in  emu/gm) 


Hence, 


and 


W 
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-216 


He  =  |  |  Mg  =  216X  3.04  X  10^  =  6.56  x  107  amp/ meter 


(or  8.27  X  105  oe) 

In  Fig.  6,  x/p  has  been  plotted  as  a  function  of  temperature  for  the  same  crystal  orientation. 
In  the  temperature  range  below  T^,  two  curves  are  shown  corresponding  to  H^  equal  to  2  and 
10  koe.  The  high  field  susceptibility,  which  is  essentially  equal  to  Xj/p,  is  constant  up  to  a 
temperature  of  80°K  and  then  increases  until  T^  is  reached.  The  low  field  susceptibility  is  equal 
to  0.765  X|/p  at  4.2°K,  climbs  with  temperature,  and  is  roughly  equal  to  Xj  /p  at  T^.  Clearly, 
in  the  low  field  case,  all  the  magnetization  cannot  be  along  the  easy  axis  parallel  to  Hj^;  other¬ 
wise,  the  susceptibility  at  4.2°K  would  be  much  smaller. 

Figure  7  shows  the  variation  of  (x/p)~  with  temperature.  The  value  of  P  calculated 
from  the  slope  of  the  curve  in  the  paramagnetic  region,  is  5.96  Bohr  magnetons.  The  negative 
temperature  intercept  is  — 160°K.  From  Eq.  (2), 

-160  =  CIW^  +  Wl2)  . 


Hence, 

W,,  +  =  -238 

11  1  - 


and  therefore, 


W 


11 


=  -22 


Curves  of  x/p  vs  H^^,  x/p  vs  temperature,  and  (x/p)  *  vs 
Figs.  8  through  10  for  the  case  of  H^^  parallel  to  a  [100]  axis, 
observed  as  for  the  previous  orientation.  The  value  of  (W^l, 


temperature  are  shown  in 
The  same  general  trends  are 
calculated  from  the  large  field 
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Fig. 4.  (x/p)  ^  vs  temperature  for  RbMnFg 
powder  specimen. 


Fig. 5.  x/p  vs  Hq£  ^or  RbMnF3  single  crystal; 
Hpc  Para^el  1 1]  axis. 


Fig.  6.  x/p  vs  temperature  for  RbMnFg  single  crystal; 
Hqc  parallel  to  [ill]  axis. 
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Fig.  7.  (x/p)  ^  vs  temperature  for  RbMnF3 
single  crystal;  parallel  to  [  1 1 1  ]  axis. 


Fig.  8.  x/p  vs  H[)c  for  RbMnF3  single  crystal; 
Hoc  parallel  to  [  100]  axis. 
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<0 


Fig.  9.  x/p  vs  temperature  for  RbMnF3  single  crystal; 
parallel  to  [  1 00]  axis. 


Fig.  10.  (x/p)  ^  vs  temperature  for  RbMnFg  single 
crystal;  parallel  to  [  1 00]  axis. 
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—  1 

susceptibility  at  4.2°K,  is  216.  The  value  of  P^,  calculated  from  the  (x/p)  vs  temperature 
plot,  is  5.8  Bohr  magnetons,  which  is  somewhat  low.  The  negative  temperature  intercept  is 
—  192°K,  giving  a  value  of  —80  for  W^. 

D.  DISCUSSION  OF  RESULTS 

The  measured  value  of  T^  was  about  94°K  for  both  the  single  crystal  and  powder  samples  of 
RbMnF^.  This  result  is  significantly  different  from  the  value  reported  by  Freiser,  et  al  2  Errors 
due  to  incorrect  location  of  the  thermocouple  junction  have  already  been  discounted.  The  liquid 
nitrogen  boiling  point  was  checked  for  calibration  purposes  by  first  using  the  ice  point  and  then 
the  boiling  point  of  liquid  nitrogen  for  the  reference  temperature.  In  both  cases,  the  calibration 
was  within  the  accuracy  (77.5  ±  1°K).  During  the  experiment,  the  sample  was  cooled  by  a  con¬ 
tinual  flow  of  gas  rising  from  the  evaporating  coolant.  Since  the  gas  had  low  heat  capacity,  it 
is  possible  that  conduction  of  heat  along  the  thermocouple  leads  raised  the  junction  temperature 
sufficiently  to  account  for  the  discrepancy. 

The  absolute  accuracy  of  the  vibrating  sample  magnetometer  ultimately  relies  on  locating 
the  sample  exactly  at  the  saddle  point.  The  same  precaution  must  be  observed  when  calibrating 
the  instrument  with  the  nickel  sphere.  It  was  estimated  that  the  absolute  accuracy  was  not  better 
than  ±3  percent.  The  relative  accuracy,  i.e.,  repeatability  of  the  measurements  once  the  sam¬ 
ple  has  been  positioned  at  the  saddle  point,  can  be  much  less  than  1  percent.  However,  con¬ 
siderable  care  must  be  taken  to  avoid  the  many  sources  of  error.  Common  problems  encountered 
were:  (1)  condensation  of  water  vapor  onto  sample  and  capillary,  (2)  mechanical  obstruction 
caused  by  the  sample  touching  the  capillary  tube,  (3)  error  signal  induced  in  the  pickup  coil  due 
to  contact  between  the  pickup  coil  and  dewar  vessel,  (4)  drift  in  the  amplifiers  due  to  fluctuating 
battery  voltages. 

There  is  excellent  agreement  between  the  calculated  and  published  value  of  P  ^  for  the 
powder  sample  and  for  the  single  crystal  with  H^^,  along  the  [111]  axis.  However,  the  calcu¬ 
lated  value  of  Hg  and  are  7  percent  lower  than  the  value  given  by  Freiser,  et  ad.,2  whereas 

the  value  of  is  considerably  lower  (—22  compared  to  —95).  It  is  important  to  note  that 

is  calculated  by  taking  the  difference  between  two  quantities  fairly  equal  in  magnitude,  so  that 
the  above  discrepancy  is  not  surprising.  The  data  calculated  from  the  measurements  with  H^^, 

along  the  [  100]  axis  are  considered  less  reliable  because  of  the  low  value  of  P  ff. 

15  en 

It  has  been  shown  that  for  the  N§el  ground  state  to  be  stable, 

whmi  i 

W^mT  mUSt  be  >"  J  ■ 

12  2 

For  the  two  pairs  of  experimental  values,  W^/W^  <  l/2.  At  4.2  °K,  =  —  M^.  Hence,  it 

follows  that  the  stability  criterion  is  satisfied.  (Note  that  the  value  of  T^  quoted  in  Ref.  1  would 
violate  the  criterion.) 

The  measurements  of  \/p  vs  for  s*ngle  crystal  do  not  indicate  any  abrupt  spin 

flopping  over  the  entire  range  of  applied  field.  The  simple  theory  outlined  in  Ch.  Ill  predicts 
that  flopping  should  occur  for  at  least  the  [111]  direction,  provided  M  and  H  are  in  the  same 
quadrant  of  the  (110)  plane,  although  misalignment  of  H^^,  would  reduce  the  effect  (see  Fig.  12). 
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CHAPTER  IH 

STATIC  EQUILIBRIUM  PROBLEM 


A.  THE  MODEL 


The  determination  of  the  static  equilibrium  of  the  magnetization  in  a  two-sublattice  anti- 
ferromagnet  is  a  four-variable  problem.  The  variables  may  be  taken  to  be  the  spherical  position 
coordinates  (0^,  i  =  1,  2,  of  each  sublattice.  The  process  of  finding  the  equilibrium  position 

involves  minimizing  the  total  free  energy  of  the  system.  For  a  cubic  antiferromagnet,  this  re¬ 
quires  the  solution  of  four  inhomogeneous  simultaneous  transcendental  equations.  The  result  is, 
in  general,  multivalued.  However,  by  restricting  H^^,  to  lie  in  the  (110)  plane,  the  problem  be¬ 
comes  much  simpler.  Despite  the  loss  in  generality,  the  restricted  solution  of  the  equilibrium 
problem  is  still  very  pertinent  to  the  analysis  of  resonance  experiments,  particularly  those  re¬ 
ported  in  Ch.  IV. 

Figure  11  illustrates  the  directions  of  H^^  and  the  sublattice  magnetizations  in  the  (110) 
plane.  The  direction  of  H^^,  makes  an  angle  ip  with  the  [100]  axis.  The  directions  of  and 
are  0  and  0  +  tt  +  A0,  respectively. 

1 3 -44 -8881  \ 


Fig.  11.  Illustration  of  M  and 
lying  in  the  (110)  plane. 


The  model  is  based  on  the  assumption  that  the  magnetization  is  locally  uniform,  i.e.,  a 
single  domain  is  considered.  A  second  assumption  is  that  Hj^,,  M^,  and  are  coplanar.  The 
statement  is  correct  for  certain  directions  of  and  is  expected  to  hold  for  a  range  of  ip  from 

symmetry  considerations.  The  fact  that  the  anisotropy  surface  is  continuous  and  has  a  stationary 
value  in  the  (110)  plane  indicates  that  the  total  energy  can  have  a  stationary  value,  which  may  or 
may  not  represent  stable  equilibrium.  However,  if  these  assumptions  lead  to  a  correct  predic¬ 
tion  of  any  of  the  observed  microwave  resonances,  they  will  be  justified. 
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B.  DERIVATION  OF  THE  EQUILIBRIUM  EQUATION 


The  total  free  energy  may  be  written  as 

K=  K  +  K  + 

an  ex  Z 

where  the  subscripts  refer  to  the  anisotropy,  exchange,  and  Zeeman  energies,  respectively. 

K 


an 


K  4  Z  Z  4  Z  Z 

-4-  (sin  ©1  sin  2<p ^  +  sin  Z©^  +  sin  ©2  sin  Z <p^  +  sin  Z©2) 


where  K  is  a  positive  constant. 


*!=*Z  =  4 


©4  =  © 


©Z  =  ©1  +  ^  +  A© 


Hence, 


K  =  —  -x*  [sin  4  ©  +  sin  2  Z©  +  sin  4(0  +  A©)  +  sin  2(Z©  +  ZA©)1 
an  4  1  1 


The  intersublattice  exchange  energy  is 

Hex=^oWl2^1  '  ^2-i^oWlllMl|2-I^oW22lM2|: 

where 


=  M  =  -M2  and  =  W22 


Assuming  that  the  lengths  of  the  vectors  M.  and  M9  remain  fixed,  the  second  two  terms  in  K 
do  not  enter  into  the  minimizing  of  the  total  free  energy  with  respect  to  the  positions  of  and 
M2  and  can  be  dropped. 

Thus, 


K 


—  fi  W  9  M.  •  M9  =  (jl  W,9M  cosA© 
no  1Z  1  Z  ro  1Z 


and 


Hz=-tXo(M1+M2)'  fi 


=  —  fjLQMH[cos  (©  —  ip)  —  cos  (©  +  A©  —  ip)] 


where  H  is  the  DC  field. 

Hence,  the  total  energy  is  given  by 


K  =  (jl o W ^  2 IVI  cos  A©  —  (l^MHfcos  (©  —  ip)  —  cos  (©  +  A©  —  ip)] 

K  4  Z  4  2 

-  4-  [sin  ©  +  sin  Z©  +  sin  (0  +  A©)  +  sin  Z(©  +  A©)] 


(3) 


By  minimizing  K  with  respect  to  A©, 

aK 


0A© 


0 


which,  after  some  manipulation,  becomes 

[Hex  —  H  cos  (©  —  ip)]  sin  A©  =  H  sin  (©  —  ip)  cos  A©  +  -g  ( — )  [2  sin  Z(©  +  A©) 

*o  VA 

+  3  sin4(©  +  A©)] 
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where 


H  =  |w,,|m  . 

ex  1  12' 

In  conformity  with  the  usual  definition  of  effective  anisotropy  field, 

= _ L  an 

an~  ^o  9M. 


and  H  =  4K/3jjl  M,  where  H  is  the  effective  anisotropy  field  along  the  [111]  axis.  By  expanding 

3-  O  3. 

the  terms  involving  A0  and  rearranging  the  equation, 

3H 


AG  = 


H  sin(0  —  ip)  +  (2  sin  20  +  3  sin  40) 


3H 

Hpx  —  H  cos  (0  —  ip)  +  — (cos  20—3  cos  40) 


We  now  examine  the  magnitude  of  the  various  quantities  in  the  above  expression.  The  largest 

value  of  H  that  will  be  considered  is  lOkoe.  Since  H  and  H  are  4.47  and  8.9  X  10  oe,  respec- 

a  ex 

tively,  there  is  clearly  negligible  error  in  the  approximation: 

AG  e  H  sin  (0  -»)  (4) 

ri 


The  term  involving  H  in  the  denominator  can  be  discarded  because,  if  H  is  large,  cos(0  —  ip)  = 
90°.  Hence,  the  term  is  always  small.  The  above  result  shows  that  even  for  the  largest  value 
of  H  to  be  considered,  A0  is  only  a  fraction  of  a  degree. 

It  is  noted  that  Eq.  (4)  does  not  contain  Ha  explicitly.  Hence,  finding  the  equilibrium  position 
first  involves  a  trade-off  between  Zeeman  and  exchange  energies,  followed  by  a  rotation  of  the 
magnetization  to  minimize  the  total  energy. 

i  2 

By  using  the  small  angle  approximations,  sinA0  =  0  and  cosA0  =  1  —  *( A0)  ,  and  neglect¬ 
ing  the  effect  of  A©  on  the  anisotropy  energy,  Eq.  (4)  enables  the  total  free  energy  equation  to  be 
written  as 


cos  (0 


At  =  -Hex  ['  -  i  (irx) 2  sin  2(0  -  *']  -H 

+  (h^-)  ^ )]  J  —  ^  H^fsin^G  +  sin^20) 


-t)  (J-)2  IsiaWli 

'  O  V  * 


(5) 


By  minimizing  K  with  respect  to  0,  (l/poM)  (8K/30)  =  0,  which  leads  to  the  result: 


4H 


3H  H 
ex  a 


Is  in  2(0  —  ip)  +  j  sin(0  —  ip)  [1  +  3  cos  2(0  —  ip)]\  =  —  sin  20(1  +  3  cos  20) 

'  ov'  I 


(6) 


The  second  term  on  the  left  side  of  Eq.  (6)  is  quite  negligible  compared  to  the  first  term,  since 
H/H  =  1/100,  except  in  the  pathological  case  of  sin  2(0  —  ip)  =  0.  Even  for  this  case,  the  re¬ 
sulting  error  in  0  is  small.  Hence,  the  final  expression  for  the  equilibrium  is 


4H 

3H 


2 

H 


ex 


a 


sin  2(0  —  ip)  =  —sin  20 ( 1  +  3  cos  20 ) 


(7) 


and  the  total  free  energy  at  equilibrium  is  given  by 


(8) 


K  _  „ 

11  M'  ex 
ro 


H 


^  sin  2(0  —  0)  — 


2H 


sin  (0  —  0)  cos(0  —  0) 


-|Ha(sin4e 


+  sin  20) 


Whether  the  equilibrium  is  stable  or  unstable  may  be  determined  by  examining  the  second 
derivative  of  the  energy  and  finding  whether  it  is  positive  or  negative.  After  differentiating  and 
simplifying, 


1  32K 

^oM  39 2 


H 

H 


cos  2(0  —  0)  — 


ir 


2h: 


cos  (0  —  0)  [2  —  9  sin  (0  —  0)] 


ex 


(cos  20  +  3  cos  40 ) 


(9) 


C.  PARTICULAR  EQUILIBRIUM  SOLUTIONS 
Case  1 :  H  =  0 

Equation  (7)  reduces  to  sin  20  (1+3  cos  20)  =  0.  Hence,  either 
sin  20  =  0  ,  i.e.,  0  =  (2n  +  1)  ^ 


or 

cos20=-|  ,  i.e.,  0  =  ±54.7°  or  ±(180°  -  54.7°)  . 

For  0  =  [(2n  +  l)7r]/2,  Eq.  (9)  shows  that  these  are  unstable  equilibria.  The  solutions  represent 
the  magnetization  along  the  cube  edges,  which  are  the  hard  directions. 

The  solutions  for  which  cos  20  =  —  ( 1  /3 )  represent  the  four  <11 1>  axes.  Equation  (9)  shows 
that  these  are  stable  equilibria. 

Case  2:  H  Along  a  <lll>_pirection  (0_=  54.7°) 

Unless  II  =  0,  as  already  discussed,  sin  20  =  0  is  no  longer  a  solution.  Two  solutions  exist 
either 

sin2(0  —  ip)  =  0  and  (1  +  3cos20)  =  0 
or 


4H‘ 


3H  H 
ex  a 


sin  2(0  —  0)  =  —  sin  20(1  +  3  cos  20) 


The  first  solution  corresponds  to  the  spins  lying  along  the  [111]  axis  and  parallel  to  H. 
Substitution  of  the  constraint  0=0=  54.7°  into  Eq.  (9)  gives  the  condition 


1 


ML  M 
^o 


ex 


II  “ 


+  2H 


H 


ex 


>  0 


2 

for  stability.  By  neglecting  the  second  term,  H  <  2HexHa. 

For  the  second  solution,  when  H  =  0,  0  =  (180°  —  54.7°);  and  when  H  =  0  =  (90°  +  ip)t  i.e., 

the  spins  are  orthogonal  to  H.  It  follows  that  125.3°  <  0  <  144.7°.  From  this  inequality,  it  can 
be  shown  that  the  second  solution  always  represents  a  stable  equilibrium. 
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Hence,  as  H  is  increased  from  zero  strength,  any  spins  initially  directed  along  the  [111] 
axis  parallel  to  H  will  remain  there  until  H  exceeds  the  value  / 2H  H  .  The  only  other  stable 
equilibrium  position  is  then  at  an  angle  lying  between  125.3°  and  144.7°,  and  the  magnetization 
will  suddenly  rotate,  or  flop,  to  the  new  position.  This  is  analogous  to  the  spin  flopping  in  a 
uniaxial  antiferromagnet.  There  is,  however,  an  important  difference.  The  spins  will  not  flop 
back  when  H  is  reduced  again,  but  will  slowly  rotate.  At  zero  field,  the  spins  will  point  in  the 
direction  0  =  125.3°,  i.e.,  along  the  other  [111]  axis  in  the  (110)  plane.  According  to  the  simple 
theory,  at  least  25  percent  of  the  magnetization  will  initially  point  in  the  latter  direction.  These 
spins  will  never  flop,  but  will  rotate  slowly  under  the  action  of  varying  field  strength  as  outlined 
above. 

Case  3:  H  Along  a  <100>  Direction  (ft  =  0) 

Solutions  to  Eq.  (7)  are 

sin  20  =  0  ,  i.e.,  0  =  (2n  +  1)  J 


4H 


3H  H 
ex  a 


=  — ( 1  +3  cos  20) 


The  solution,  sin  20  =  0,  is  stable  if  H  >  3H  H  /2. 

'  ex  a 

2 

The  second  solution  becomes  imaginary  if  H  >  3H  H  /2, .  but  for  lower  values  of  H,  it  is 
the  only  solution  and  is  stable. 

Case  4:  H  Along  a  <  1 1 0>  Direction  (rl>  =  90°) 


Equation  (7)  becomes 
t2 


4H 


3H  H 
ex  a 


sin  20  =  sin  20(1  +  3  cos  20) 


The  two  solutions  are  either 
sin  20  =  0 


i.e.,  0  =  mr 


or 


_ 4H 

3H 


2 

H 


ex 


a 


=  1+3  cos  20 


The  first  solution  is  stable  if  H2  >  3H  H  . 

ex  a  ^ 

The  second  solution  does  not  exist  if  H  >  because  cos  20  is  greater  than  unity. 

For  the  second  solution  to  be  stable. 


I  H2  3"»'l  M 

Ih - 5-  J  cos  20  - 

'ex  ' 


9H 


cos  40  >  0 


The  inadequacy  of  the  model  may  be  noted  by  considering  the  effect  of  a  strong  field  applied 
along  a  </110^>  direction.  The  above  result  states  that  the  magnetization  will  be  perpendicular  to 
H  and  directed  along  the  [100]  axis.  However,  the  [100]  axis  is  a  hard  axis,  and  the  spins  are 
free  to  rotate  out  of  the  (110)  plane  that  contains  H,  while  maintaining  the  orthogonality  condition 
with  respect  to  H,  until  they  lie  along  the  [110]  direction.  Clearly,  this  is  a  lower  energy  state. 
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Case  5:  H  in  an  Arbitrary  Direction 


For  the  general  case,  the  transcendental  Eq.  (7)  must  be  solved,  and  the  solution  must  then 
be  checked  for  stability.  This  was  done  for  five  values  of  ip.  The  results  are  shown  in  Fig.  12. 
The  dotted  portions  of  the  curves  for  ip  =  10°  and  ip  =  45°  represent  unstable  regions. 


(2Hei  H0  ),/2 


Fig.  12.  Equilibrium  position  of  sublattice  magnetizations  as  a  function  of  H 


DC 


D.  CONCLUSIONS 

The  model  enables  a  simple  expression  for  the  equilibrium  position  of  the  magnetization  to 
be  found.  It  appears  to  give  reasonable  solutions  for  H  parallel  to  the  [100]  and  [111]  axes,  but 
leads  to  a  contradictory  result  for  a  large  field  applied  along  a  [110]  axis. 

The  magnetic  field  acts  differently  on  the  magnetization  depending  on  whether  M  and  H  are 
in  the  same  quadrant  of  the  { 1 1 0}  plane  or  in  different  quadrants.  For  the  latter  case,  as  H  is 
increased,  the  magnetization  rotates  away  from  H.  For  very  large  fields,  M  and  H  are  perpen 
dicular.  If  M  and  H  are  initially  in  the  same  quadrant,  with  H  along  the  [111]  axis,  M  will 
keep  its  orientation  along  the  same  [111]  axis  until  H  exceeds  a  certain  value,  equal  to  J 
when  M  will  abruptly  flop  into  the  adjacent  quadrant.  If  M  and  H  are  initially  in  the  same 
quadrant  and  H  is  not  along  the  [111]  axis,  M  gradually  rotates  away  from  H  as  it  is  increased. 
After  a  certain  rotation,  depending  on  the  magnitude  and  orientation  of  H,  M  will  flop  into  the 
adjacent  quadrant  and  remain  there. 

The  simple  model  predicts  that  under  the  action  of  a  DC  magnetic  field,  a  single  crystal  of 
RbMnF^  can  act  like  a  two-state  device.  If  a  strong  field  is  applied  to  one  quadrant  in  the  (110) 
plane,  it  will  cause  any  magnetization  lying  along  an  axis  within  that  quadrant  to  rotate  into  the 
adjacent  quadrant  in  the  plane  and  remain  there  when  the  field  is  removed.  If  the  field  is  then 
applied  within  the  quadrant  to  which  M  rotated,  it  will  cause  M  to  rotate  back  to  its  original 
quadrant  and  remain  there. 
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CHAPTER  IV 

ANTIFERROMAGNETIC  RESONANCE 


A.  RESONANCE  ANALYSIS 

The  equation  of  motion  for  the  two-sublattice  antiferromagnet  is 


i,j  =  1,2  and  i  ^  j 


The  capital  letters  represent  the  steady  state  values,  the  deltas  refer  to  increments  about  the 
equilibrium  position,  and  T  represents  the  torque  resulting  from  the  anisotropy  energy.  De¬ 
magnetizing  effects  are  neglected. 

+  H^,  w^ere  nuc-^ear  hyperfine  field  and  acts  in  the  direction  of  IVL. 

The  nuclear  hyperfine  interaction  is  important  for  RbMnF^,  since  at  low  temperatures,  it  is 
comparable  in  magnitude  [(-  9.46/T°K)  oe]  to  the  effective  crystalline  anisotropy  field.  The  terms 
that  contribute  to  the  static  equilibrium  condition  may  be  separated  from  the  dynamic  terms: 


M. 

— -  =  0  =  u.  M.  X  (H  .  +  W..M.  +  W..M.)  +  T. 
y  ro  i  oi  n  l  ij  j  l 

Therefore, 


<5M. 

- -  =  |j.  [<5M.  X  (H  .  +  W..M.)  +  M.  X  W..6M.]  +  ST. 

Y  o'  i  oi  i]  j  l  lj  J  i 

The  anisotropy  energy  for  the  two-sublattice  cubic  system  can  be  expressed  as 


(10) 


and 


H  = - (M2  M2  +  M2  M2  +  M2  M2  +  M2  M2  -  +  M2  M?  +  M?  M2  ) 

an  ^j4  lx  ly  lx  lz  ly  lz  Zx  Zy  Zx  Zz  Zy  Zz 


H  .  =__L 

aix,y  ^o  8M. 


x,y 

neglecting  intersublattice  anisotropy.  The  anisotropy  torques  are  found  by  taking  the  vector 

product  M.  X  H  ..  Assuming  that  the  axes  of  the  Cartesian  coordinate  system  are  chosen  to  coin- 
i  a  i 

cide  with  the  cube  edges. 


ZK"  7  7 

T.  =  M.  H  .  -M.H.  =  M.  M.  (M  -  M  ) 


ix  ly  aiz 


iz  aiy 


M 


ly  iz'  ly 


Hence, 


T.  =  M.  H  .  -  M.  H.  =  -^  M.  M.  (M2  -  M2  ) 

iy  iz  aix  ix  iz  ^4  ix  iz  iz  ix 

T.  =  M.  H  .  -M.  H.  =  M.  M.  (M?  -  M2  ) 

iz  IX  aiy  iy  IX  M4  ix  iy  ix  iy 


T.  +  ST.  =  (M.  +  6M.  )  (M.  +  5M.  )  [(M.  +  6M.  )2  -  (M.  +  6M.  )21 

ix  IX  tv/t4  iy  IV  iz  iz  lx  IV  IV  iz  iz  ] 


M 


iy  iz 


iy 


iy 


zi 


It  follows  that 


Similarly, 


and 


6T.  =  ^  [6M.  (3M.  M?  -  M?  )  -  6M.  (3M?  M.  -  M3  )] 

ix  m4  1  ly  iz  ly  iz  iz  iz  ly  ly  J 

6T.  =  [6M.  (3M?  M.  —  M?  )  —  6M.  (3M?  M.  —  M?  )] 

iy  ^4  1  iz  iz  ix  ix  ix  ix  iz  iz  1 

6T.  =  ^  [6M.  (3M?  M.  -  M?  )  -  5M.  (3M2  M.  -  M3  )] 

iz  ^4  1  ix  ix  iy  iy  iy  iy  ix  ix  J 


Equation  (10)  is  now  linearized  and  written  in  terms  of  Cartesian  components.  Six  simul¬ 
taneous  equations  are  obtained;  the  variables  are  the  six  differential  increments  <5M.  ,  6M.  , 

ix  iy 

In  matrix  notation,  the  set  of  equations  may  be  written  as 


|6M]  =  [A]  [6M]  +  [B]  [6M] 


(11) 


Matrix  [A]  contains  the  coefficients  of  6M  which  are  invariant  to  a  rotation  of  the  coordinates; 
[B]  contains  the  anisotropic  terms  which  are  affected  by  a  rotation. 

The  coordinates  are  now  rotated  through  an  angle  0  in  the  (110)  plane  so  that  the  new  z-axis 
and  the  equilibrium  direction  of  the  magnetization  coincide.  The  new  coordinates  will  be  re¬ 
ferred  to  as  (x,y, z).  Thus,  M  =  M,  and  M  =  M  =0. 

J  z  x  y 

If  [T]  is  the  transformation  matrix  from  the  new  to  the  old  axes  [6M]  =  [T]  [6m],  where  [6m] 
represents  the  dynamic  components  of  the  magnetization  in  the  new  coordinate  system.  Hence, 
from  Eq.  (11), 

3Y  =  [A]  l6ml  +  [T]"1  IB]  m  l6m] 

where  the  components  in  [A]  are  written  in  terms  of  the  new  coordinates.  The  equation 

[6m]=[a][«m]  +  [b|[6ml  (12) 

defines  the  new  matrices  [a]  and  [b ]. 

The  derivation  of  the  transformation  matrix  is  given  in  the  Appendix.  For  the  two-sublattice 
system,  it  becomes  a  6  X  6  matrix: 
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m  = 


.  cos  0  +  1  '  ^  cos  6  —  1  j  sine  Q 

2  2  41 


^  cos  9  —  1  j  (  cos_0_+_l  j  sin  9  Q 


—  sin  9 

42 


—  sin  9 

4z 


4z 


cos  e 


(  COS0jt_l  ) 


(  cose  -  1 ) 


—  sin  9 


42 


(  cose  -  1 ) 


( cosej_i ) 


—  sin  9 


41 


sin  9 

42 

sin9 

42 


cos  9 


and 


[T]'1  =  [Tl1 


By  invoking  the  constraint  that  the  total  magnetization  must  be  conserved,  [a]  and  [b]  may 
be  reduced  from  6  X  6  to  4  X  4  matrices. 


m 2  =  M2  =  (m  .  +  6m,  )2  +  (m,  +  6m,  )2  +  (m  +  6m,  )2 

1  lx  lx  ly  ly  lz  lz 


Putting 


2  ,  2  2  _  ,.2 

m ,  +  m.  +  m .  ~  M 

lx  ly  lz 


and  neglecting  the  second  order  terms 


m.  6m,  +  m.  6m,  +  m ,  6m.  =  0 

lx  lx  ly  ly  lz  lz 


Since  mlx  =  miy  =  an<^  miz  ^  ^  follows  that  <5mlz  -  0.  Similarly,  ^m2z  ~  Hence, 

the  rows  and  columns  in  [a]  and  [b]  associated  with  6mlz  and  Sm^  may  be  struck  out. 

After  considerable  algebraic  work. 


[al 


(H  ~  +Hm  +  H  ) 
DCz  N  ex 


‘^DCz^N^W 


H 


-H 


-H 


H 


<HDCz-HN-»eX' 


-|HDC2-HN-Hcx' 
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[b] 


bll  bl 2 


~bl  2  ~bll 


bll  b12 


0  °  bl  2  bU 


and 

,  3K  .  2„  /  2„  sin20\ 

bn  =  ]Tm  sin  e  lcos  e - —) 

u  2K  /,  .  2  2„  ,  sin40\ 

b12  =  TT"m  \3  Sln  e  cos  6  “  cos  0  +  4  ) 

Equation  (12)  is  rewritten  as 

[a  +  b  —  ^  ]  [6m]  =0  or  [c]  [6m]  =  0 

The  determinant  of  [c]  must  be  equal  to  zero  for  the  set  of  linear  equations  to  be  simultane¬ 
ously  true.  By  multiplying  out  the  determinant,  a  biquadratic  equation  results  in  u>/y: 

D(  )4  +  E(— )2  +  F  =  0 
Wo  wo 

where  D,  E,  and  F  are  functions  of  the  elements  of  [a]  and  [b]. 

Hence, 


)2 


-E 


—  4DF 


After  the  appropriate  substitutions  for  D,  E,  and  F  and  considerable  algebra,  the  following 
expression  for  the  normal  mode  resonant  frequencies  is  obtained. 

=  HDCz  +  2IWHN  +  b12J  +  (HN  +  b12)  _bll 

1  ro 


2HDCz  ./2H„|H  N  *  b12>  *  |HN  *  b12|2  + 


(13) 


The  equation  yields  only  two  independent  frequencies.  It  is  interesting  to  note  that  co  is 
independent  of  and  It  is  only  a  function  of  the  component  of  along  the  equi¬ 

librium  position  of  M. 


B. 


NORMAL  MODE  FREQUENCIES  FOR  HDC  PARALLEL  TO  A  [111]  AXIS 

For  this  important  case,  b^  =0,  and  =  Ha-  ^'cluation  (13)  becomes  a  perfect  square: 


25 


or 


«#'2  ■  <HDCz  *  '2Hex(HN  +  "a1  *  («N  *  h/^ 

1  ro 

^  =  Hdc  cos  (0  -  *)  *  J  (Hn  +  Ha)  (2Hex  +  Hn  +  Ha)  .  (14) 

Now,  if  <  J  2H^xHa,  one  of  the  permitted  equilibrium  conditions  is  0  =  ip  =  54.7°.  Then 

Eq.  (14)  is  equivalent  to  the  well-known  result  of  Keffer  and  Kittel.^  However,  there  is  another 
pair  of  frequencies  that  corresponds  to  the  magnetization  which  was  initially  directed  along  the 
second  [111]  axis  (not  parallel  to  Hj^).  The  resonant  frequencies  for  both  pairs  of  modes  have 
been  plotted  as  a  function  of  Hj^  in  Fig.  13.  The  low  frequency  branch  should  approach  zero 
frequency  at  the  spin  flop  field.  The  simple  theory  has  neglected  the  nuclear  magnetic  reso¬ 
nance  (NMR)  modes  which  exist  at  much  lower  frequencies  than  the  electronic  modes.  At  low- 
frequencies,  the  electronic  and  nuclear  spin  systems  will  couple  in  such  a  way  that  one  branch 
of  the  coupled  mode  spectrum  will  go  to  zero  frequency  at  the  spin  flopping  field. 


C.  RESONANCE  EXPERIMENTS 

Resonance  experiments  were  performed  at  X-band  and  at  4.2°K  with  the  experimental  ar¬ 
rangement  illustrated  in  the  block  diagram  of  Fig.  14.  The  single  crystal  RbMnF^  specimen 
was  a  4-mm  cube  which  was  cut  so  that  four  of  its  faces  were  (110)  planes  and  two  were  (100) 
planes.  The  crystal  was  mounted  in  a  TE^qx  cavity  at  a  point  of  maximum  RF  magnetic  field 
and  could  be  rotated  about  a  [110]  axis.  Three  different  cavities  were  used  so  that  experiments 
could  be  performed  at  five  frequencies  distributed  over  the  X-band  spectrum.  The  cavity  and 
connecting  waveguide  were  immersed  in  liquid  helium.  A  glass  double  dewar  system  enabled 
the  sample  to  be  maintained  at  4.2°K  for  many  hours  at  a  time.  The  dewar  system  was  sus¬ 
pended  between  the  pole  pieces  of  a  12-inch  Varian  electromagnet.  The  waveguide  and  cavity 
were  positioned  so  that  the  RF  magnetic  field  at  the  crystal  and  Hj-^  were  orthogonal.  Rotation 
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Fig.  14.  Block  diagram  of  X-band  AFMR  experimental  arrangement. 

of  the  crystal,  which  permitted  the  orientation  of  H.^  with  respect  to  the  cubic  lattice  to  be 
varied  during  the  experiment,  was  obtained  through  a  gear  train  mounted  on  the  end  of  the  cavity. 

A  shaft  rotation  was  translated  to  the  sample  mount  inside  the  cavity  by  a  thin  dielectric  rod.  It 
was  necessary  to  evacuate  the  waveguide  and  dewar  vessel,  and  refill  it  with  gaseous  helium  be¬ 
fore  transferring  liquid  helium,  in  order  to  prevent  the  gear  train  from  sticking. 

Each  of  the  cavities  was  of  the  reflection  type  and,  by  appropriate  choice  of  iris,  was  ar¬ 
ranged  to  be  overcoupled  at  4.2°K.  Resonance  was  detected  by  the  decrease  in  reflected  power 
from  the  cavity.  The  power  incident  on  the  sample  was  about  lOmw.  The  resonance  linewidth 
did  not  change  when  the  incident  power  was  reduced  to  1  mw,  indicating  that  saturation  effects 
were  absent.  Reflected  power  vs  was  plotted  for  the  field  range  0  to  lOkoe  for  a  large 

number  of  crystal  orientations,  but  with  restricted  to  the  (110)  plane.  It  was  found  that 

some  of  the  observed  resonances  were  very  sensitive  functions  of  angle;  in  order  to  keep  track 
of  their  behavior,  the  angle  increment  had  to  be  about  five  degrees. 

D.  EXPERIMENTAL  RESULTS 

The  DC  magnetic  field  required  to  excite  resonance  has  been  plotted  as  a  function  of  field 
direction  for  the  five  X-band  frequencies.  The  results  are  shown  in  Fig.  15(a-e).  It  is  evident 
that  in  all  cases,  more  than  one  resonance  was  observed.  For  the  two  lowest  frequencies,  there 
are  many  resonances.  Most  of  the  resonances  can  be  tracked  over  a  wide  range  of  angle.  By 
comparing  the  shapes  of  the  curves,  it  is  possible  to  correlate  some  of  the  resonances  as  a 
function  of  frequency.  In  this  manner,  it  has  been  possible  to  construct  plots  of  resonant  fre¬ 
quency  vs  Hqq  for  two  of  the  symmetry  axes  of  the  crystal  —  the  [111]  and  [100]  axes.  The  data 
for  the  [111]  axis  have  been  superimposed  on  the  theoretical  curves  of  Fig.  13  for  comparison. 

It  is  interesting  to  note  that  there  are  no  experimental  points  corresponding  to  the  case  of  M 
parallel  to  H.  Secondly,  the  experimental  points  are  in  reasonable  agreement  with  the  resonance 

theory,  but  correspond  to  the  spins  lying  near  to  the  other  [111]  axis  in  the  (110)  plane.  The 

2 

theoretical  curve  for  the  field  dependent  mode  reported  by  Freiser,  et  al.,  has  also  been  ap¬ 
pended.  It  does  not  appear  possible  to  excite  this  mode  below  about  12Gcps  for  the  [111]  axis. 

The  graph  of  resonant  frequency  vs  for  parallel  to  the  [100]  axis  is  shown  in  Fig.  16. 

Many  more  resonances  were  observed  for  this  case  than  for  the  [111]  case.  Again,  five  of  the 
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Fig.  15(a-e). 


DC  magnetic  field  required  for  resonance  vs  angle. 
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Fig.  15.  Continued 


29 


(koe) 


FREQUENCY  (Gcps) 


(e) 


Fig.  15.  Continued. 


Fig.  16.  Observed  resonances  for  parallel  to  [  1 00]  axis. 
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resonances  seem  to  agree  fairly  well  with  the  resonance  theory.  In  this  case,  the  theory  pre¬ 
dicts  only  one  pair  of  mode  branches.  By  studying  the  curves  in  Fig.  15(a-e),  it  has  been  possi¬ 
ble  to  group  most  of  the  remaining  points  into  two  families  which  seem  to  belong  to  different 
modes  of  resonance.  One  group  falls  very  close  to  the  theoretical  plot  of  Freiser's  field  depend¬ 
ent  mode.  The  second  group  seems  to  fall  on  a  line  parallel  to  the  first  group. 

The  intercept  of  the  curve  of  AFMR  frequency  vs  Hnr  on  the  frequency  axis  is  nearly  equal 
1  /  2  U{^ 

to  ypo[2Hex(H^  +  Ha)]  •  Hence,  extrapolation  of  the  experimental  data  provides  an  estimate 

of  this  quantity.  Unfortunately,  there  are  only  five  data  points,  but  an  extrapolation  indicates 

that  the  intercept  is  at  9.3Gcps.  This  is  to  be  compared  with  a  frequency  of  9.65  Gcps  which  was 

2 

calculated  by  using  the  data  of  Freiser,  et  al.  The  discrepancy  can  be  accounted  for  by  a  7- 

percent  difference  in  the  value  of  H  (  or  a  difference  in  H  ).  Note  that  this  is  consistent  with 

ex  a 

a  similar  discrepancy  noted  in  the  susceptibility  data. 

The  shape  of  the  plots  of  resonance  field  vs  direction  of  applied  field  gave  an  excellent  check 
on  the  accuracy  of  the  orientation  of  the  crystal,  since  the  curves  should  have  mirror  symmetry 
about  the  zero  and  90°  points.  The  curves  in  Fig.  15(a-e)  have  been  replotted  with  the  angular 
corrections,  which  did  not  exceed  3°  in  any  of  the  experiments. 
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CHAPTER  V 

EXCITATION  OF  SPIN  WAVES  BY  PARALLEL  PUMPING 


A.  MEASUREMENT  OF  SPIN  WAVE  LINE  WIDTH  BY  PARALLEL  PUMPING 

Spin  wave  instability  and  premature  saturation  of  the  AFMR  have  been  analyzed  by  Heeger 
and  Pincus.  According  to  their  theory,  the  instability  threshold  is  given  by 


/W«AHIr  X1/2 


n*  c>  kv 

h  .,  =  4AH  [ — ^ - -) 

°\  / 


crit 


where  AH  is  the  uniform  mode  linewidth,  AH,  is  the  spin  wave  linewidth,  and  a)  is  the  AFMR 
o  k  o 

frequency. 

9 

Saturation  of  AFMR  in  KMnF^  has  been  observed  by  Heeger  who  found  that  the  onset  of 
saturation  occurred  at  a  very  low  threshold  —  on  the  order  of  1  mw.  He  inferred  that  the  spin 

_3 

wave  linewidth  was  on  the  order  of  5  X  10  oe,  compared  with  40  oe  for  the  uniform  mode. 

Naiman  and  Lawrence1  ^  have  observed  a  very  long  relaxation  time  for  spin  flopping  in  the  anti- 
ferromagnet  CuCl^  H^O  and  have  compared  their  results  to  Heeger' s.  The  question  arises 
whether  similar  saturation  effects  occur  for  RbMnF^.  Thus,  it  was  proposed  that  the  spin  wave 
linewidth  for  RbMnF^  be  estimated  by  measuring  the  threshold  RF  magnetic  field  for  longitudi¬ 
nally  pumped  spin  waves. 

17  18 

Morgenthaler  and  Schlomann  first  proposed  that  spin  waves  at  half  the  pumping  frequency 

could  be  parametrically  excited  by  the  application  of  an  RF  driving  field  parallel  to  the  DC  mag- 

19 

netizing  field.  Morgenthaler  has  extended  his  analysis  to  include  two-sublattice  systems  and, 
in  particular,  has  derived  an  expression  for  the  instability  threshold.  1 

First  order  nonlinear  coupling  processes  are  assumed  to  exist  so  that 


=  Iwj  I  +  |w2| 


and 


kl  =  ~k2 


The  terms  u)  ^  and  o)2  are  the  two  magnon  frequencies,  and  k2  are  the  respective  wave  vectors, 
and  cjp  and  k^  are  the  corresponding  pump  parameters.  (Higher  order  processes  can  exist,  but 
these  generally  have  higher  thresholds  than  first  order  processes.)  Thus,  the  net  magnetization 
of  each  lattice  may  be  represented  by  a  steady  state  component  plus  a  small  component  due  to  a 
standing  spin  wave.  The  equation  of  motion  has  additional  terms  compared  to  the  uniform  pre¬ 
cession  case  (k.  =  0)  which  arise  because  of  the  nonzero  divergence  of  M.  Also  included  is  a 
term  due  to  the  RF  driving  field,  and  a  damping  term  of  the  Landau- Lifshitz  form.  The  analysis, 
based  on  a  perturbation  technique,  shows  that  the  pump  acts  to  overcome  the  loss  in  the  system. 
To  a  first  approximation,  it  does  not  shift  the  normal  mode  frequencies,  which  are  found  by 
setting  the  pump  amplitude  and  loss  equal  to  zero  in  the  equation  of  motion. 

Determination  of  the  critical  field  for  the  onset  of  instability  involves  finding  a  solution  to 
the  equation  of  motion  for  which  the  pump  provides  just  sufficient  power  to  overcome  the  loss 
terms.  The  two  possible  frequency  relationships,  corresponding  to  degenerate  and  nondegener¬ 
ate  magnon  instabilities,  are 


«1  =  w2 


ki 


and 


Cd  .  =  0) ,  . 
l  ki 


Since  the  nondegenerate  case  can  be  expected  to  have  a  higher  threshold,  we  shall  consider  only 
degenerate  magnon  instabilities.  For  the  antiferromagnet,  it  is  found  that 
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4AHkwk  /" a  +  2oJex  +  (X11  +  X12*  k  +  wMsln  ^ 


crit 


wMSin  ^ 


wa  +  (X1 1  —  X12^  k  WM 


(15) 


and  the  normal  mode  frequencies  are 

2  221 /2 
wk  =  {(a,a  +  (X11  _X12)  k  "Ml  [a,a  +  2“ex  +  (X11  +Xl2)k  wM+CJMsin  *1>  '  ±UH 


(16) 


=  a  +  (Xll“X12)  *  WH 

where  2AH^  is  the  spin  wave  linewidth, 

CJ 

wk  =  2  P 

co  =  y\i  W, 
ex  rro  12 


WH  =“^oHDCz 

"m  =-w0m 


co  =  — vul  H 
a  r  ro  a 


ip  is  the  angle  between  the  propagation  vector  k  and  the  z-axis,  and  X^,  are  the  coefficients 

of  the  next  higher  order  terms  in  the  molecular  field  expansion.  The  analysis  assumes  that  the 
steady  state  component  of  the  magnetization  is  parallel  to  the  easy  axis  and  that  the  parallel 
susceptibility  is  zero.  Since  appears  in  the  numerator  of  Eq.  (15),  it  is  evident  that  the  mode 
with  the  lower  frequency  has  the  lower  threshold.  For  the  lower  frequency,  the  threshold  may 
be  expressed  as  the  difference  of  two  quantities: 


4AH 


crit 


w  M  sin 


Z  lwa  +  2wex  +  ^Xll  +  *T2^  k  WM  +  wMsin  ^ 


_ _H _ 

ex  I  ,  77  x  ,2 

L  nK  +  (XU  -XiZ}  k  ^MJ 


(18) 


The  first  quantity  increases  with  k.  The  second  quantity  decreases  with  k,  since  its  denomi¬ 
nator  is  a  more  sensitive  function  of  k  than  the  numerator.  Hence,  it  follows  that  the  magnons 
for  which  k  =  0  have  the  lowest  threshold.  Further,  Eq.  (18)  predicts  that  for  k  =  0  magnons, 
h  ^  is  zero  when  Hj^  =  J  2H^H~,  which  is  the  condition  for  spin  flopping  derived  in  Ch.  III. 
However,  by  using  the  more  exact  expression  for  given  by  Eq.  (16),  it  is  found  that  hcr^ 
approaches  a  finite  limit  equal  to  2AH^  as  Hj^  approaches  (but  is  less  than)  (yJ2HaHex> 

In  summary,  the  threshold  field  is  lowest  for  k  =  0  magnons  propagating  in  the  x-y  plane 
(ip  =  90°).  It  can  be  reduced  to  a  small  value  by  adjusting  the  DC  magnetic  field  to  be  almost 
equal  to,  but  less  than,  the  flopping  field.  Because  of  the  dependence  of  on  Hj^,  this  implies 
that  the  frequency  of  the  pump  is  constrained  to  be  correspondingly  low.  Measurement  of  the 
threshold  field  provides  a  direct  method  of  determining  AH^,  which  enters  into  the  threshold 
equation,  provided  that  the  other  parameters  are  known. 

In  order  to  predict  the  required  RF  power  and  frequency  for  such  an  experiment  with  RbMnF^, 

the  values  H  and  H  were  taken  from  Ref.  2.  The  exchange  constants  W.  .  and  W.  ~  were 
a  ex  b  1112 
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calculated  from  the  data  to  be  —95  and  —233,  respectively.  The  exchange  constants  A^  and  A^ 
were  obtained,  starting  from  the  expression  for  the  exchange  energy: 


X 

ex 


-2J..S.  •  S. 
ij  1  J 


where  is  the  exchange  integral.  By  expanding  and  into  Taylor  series  and  assuming  that 
only  nearest-neighbor  interactions  are  important,  it  may  be  shown  that 


A 


11 


W11L 


12 


W 


where  L  is  the  lattice  constant. 

The  spin  wave  manifold  has  been  plotted  in  Fig.  17  for  a  fixed  value  of  Hj-^.  The  shape  is 
similar  to  the  spin  wave  manifold  for  a  ferrimagnetic  material,  although  the  curves  are  nearly 
degenerate  for  small  values  of  k.  Also,  an  increase  in  shifts  the  curves  parallel  to  the 

axis  toward  the  origin  rather  than  away  from  the  origin  as  in  the  ferrimagnetic  case.  This  is 
illustrated  in  Fig.  18,  where  a  family  of  dispersion  curves  for  different  values  of  has  been 

plotted.  Only  the  curves  for  si nif)  =  1  are  shown;  AH^  was  assumed  to  be  1  oe. 

The  propagation  factor  was  computed  for  various  values  of  for  a  given  applied  DC  field, 
by  using  Eq.  (17).  Corresponding  values  of  h  ^  were  then  calculated  by  using  Eq.  (15)  and  are 
plotted  in  Fig.  19.  It  is  seen  that  hcji.^  increases  very  rapidly  with  frequency  and  the  applied 
power,  which  is  proportional  to  (h  ^)^,  becomes  prohibitively  large.  Two  curves  are  shown 
corresponding  to  two  values  of  AH^  (1  and  0.01  oe).  Only  the  curves  for  sin  iJj  =  1  have  been 
plotted  since  these  have  the  lowest  threshold.  The  value  of  was  chosen  to  be  within  5  oe  of 

spin  flop  in  order  to  limit  h  ^  to  reasonably  attainable  values  in  the  VHF  region  of  the 
spectrum.  This  choice  could  be  unreasonable  on  two  counts:  (1)  the  requirement  of  extreme 


Fig.  17.  Spin  wave  frequency  vs  spin  wave  number. 
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Fig.  18.  Spin  wave  frequency  vs  k. 


Fig.  19.  h  vs  pump  frequency. 
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stability  in  the  applied  DC  field,  (2)  dislocations  in  the  crystal  could  cause  Ha  and  to  be 
nonuniform,  resulting  in  an  ill-defined  value  of  flopping  field. 

The  choice  of  the  operating  temperature  at  which  to  perform  parallel  pumping  experiments 
is  quite  important  because  the  spin  wave  linewidth  exchange  and  anisotropy  fields  all  vary  with 
temperature  and  affect  h^.^  correspondingly.  The  most  important  variable  is  the  anisotropy 
which,  according  to  Ref.  2,  is  one-hundredth  of  its  0°K  value  at  77°K.  Hence,  considering  only 
the  variation  of  H^,  hcr^  at  liquid  nitrogen  temperature  will  be  ten  times  its  value  at  liquid 
helium  temperature. 

B.  MEASUREMENT  TECHNIQUE 

From  the  foregoing  analysis,  it  was  concluded  that  by  setting  the  DC  field  to  within  a  few 

oersteds  of  the  flopping  field,  the  instability  threshold  would  be  attained  at  VHF  frequencies 

with  modest  power  levels  if  AH^  were  on  the  order  of  0.1  oe.  Accordingly,  a  parallel  pumping 

experiment  was  designed  to  operate  at  a  frequency  of  30Mcps.  This  experiment  was  a  lumped 

circuit  version  of  the  usual  microwave  parallel  pumping  experiment,  its  principle  of  operation 

being  the  detection  of  a  sudden  increase  in  xM  when  the  instability  threshold  is  exceeded.  The 

sample  was  placed  on  the  axis  of  a  solenoid  which  was  part  of  an  RF  bridge  circuit.  An  increase 

in  the  loss  of  the  sample  would  then  cause  an  unbalance  of  the  bridge.  This  type  of  measurement 

22-24 

is  similar  to  the  measurement  of  NMR.  Two  NMR  bridge  circuits  were  considered.  The 

22 

first,  used  by  Bloembergen,  et  al.,  utilizes  a  circuit  in  which  the  responses  of  two  identical 
shunt-tuned  coils,  one  containing  the  sample,  are  balanced  against  each  other.  The  phases  and 
amplitudes  of  eachhalf  of  the  bridge  are  adjusted  for  a  null  at  the  summing  point.  The  three 
disadvantages  of  the  scheme  are:  (1)  to  avoid  excessive  drift,  the  reference  channel  should  have 
the  same  environment  as  the  sample  channel,  (2)  one  of  the  channels  must  contain  a  x/2  cable 
to  invert  the  phase,  or  a  transformer  must  be  used,  (3)  since  the  amplitude  and  phase  adjust¬ 
ments  are  not  independent,  the  adjustment  procedure  is  very  tedious. 

25 

In  the  experiment,  the  detection  scheme  used  the  bridged-T  circuit  illustrated  in  Fig.  20. 
The  sample  coil  L  is  situated  in  one  arm  of  the  bridge  which  is  nonresonant.  The  bridge  is 
balanced  by  two  independent  adjustments:  The  phase  balance  is  obtained  by  varying  the  capacitor 
C;  the  amplitude  balance  is  determined  by  the  variable  capacitor  C’.  The  conditions  for  balance 
are  independent  of  the  source  and  load  impedances. 


Fig.  20.  Bridged-T  circuit. 
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Fig.  21 .  Block  diagram  of  parallel  pumping  experiment 


Although  the  bridged-T  circuit  is  inherently  a  narrow  band  device,  this  was  not  important, 
apart  from  stability  considerations,  since  the  experiment  was  performed  at  a  single  frequency. 
The  practical  limit  of  sensitivity  achievable,  which  is  a  function  of  the  depth  of  the  null,  is 
limited  to  about  —  60  db  because  of  the  difficulty  of  maintaining  the  balance. 

The  conditions  for  balance  are 

R  =  w  C1C2  J1  +  C^]  Rl  (19) 


and 


1 

w2LC 


(20) 


where  R  is  the  equivalent  shunt  loss  of  the  coil,  and  C  is  the  shunt-tuning  capacitor.  The  Q 
of  the  inductor  is  given  by 


A  block  diagram  of  the  parallel  pumping  experiment  is  shown  in  Fig.  21.  The  sample  coil 
was  immersed  in  liquid  helium  and  located  between  the  pole  pieces  of  a  6-inch  Varian  electro¬ 
magnet.  Several  different  coil  designs  were  used  in  an  effort  to  improve  the  filling  factor,  but 
each  had  an  inductance  of  about  1  \xh  and  required  about  20  to  25  pf  of  shunt  capacitance  for  bal¬ 
ancing  purposes. 

The  bridged-T  circuit  was  located  immediately  above  the  dewar  vessel  in  order  to  minimize 
the  length  of  the  connection  to  the  coil.  A  special  low  capacitance  coaxial  transmission  line  was 
used  for  the  connection.  The  center  conductor  was  a  fine  wire  supported  by  polystyrene  foam; 
the  outer  conductor  was  stainless  steel,  copper  plated  on  the  inside  to  a  thickness  of  a  skin 
depth.  When  cooled,  the  Q  of  the  coaxial  line  plus  coil  was  320.  The  signal,  which  could  be 
square  wave  or  pulse  modulated,  was  obtained  from  a  stable  generator  at  the  milliwatt  level  and 
was  amplified  to  a  level  of  about  12  watts  by  a  class  C  tuned  power  amplifier.  The  out-of¬ 
balance  signal  from  the  bridge  circuit  was  amplified  by  two  IF  strips  each  having  30 db  of  gain. 
The  output  was  displayed  directly  on  an  oscilloscope. 

Because  of  the  sensitivity  of  the  equipment,  an  increase  in  loss  of  the  sample,  equivalent 
to  a  sample  Q  of  about  1500,  could  be  detected.  The  sensitivity  to  reactive  changes  was  about 
one  part  in  2500  (0.01  pf  change  in  C').  This  corresponded  to  a  pn  of  0.01  or  a  p'  of  5  x  10  ^  for 
0.25  gram  of  RbMnF^  with  the  coil  wound  on  the  sample. 

The  orientation  of  the  sample  was  chosen  so  that  a  [111]  axis  was  parallel  to  the  axis  of  the 
solenoid  which  was  also  parallel  to  H.^. 

C.  EXPERIMENTAL  RESULTS 


With  the  coil  and  sample  cooled  to  4.2°K  and  the  signal  power  set  for  maximum  output,  the 
DC  magnetic  field  was  increased  exceedingly  slowly  over  the  range  0  to  10koe.  A  slight  reactive 
unbalance  was  observed  which  was  attributed  partly  to  the  increase  of  x'  with  and  partly  to 

a  spurious  effect  which  was  observed  even  in  the  absence  of  the  sample.  There  was  no  unbalance 
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which  could  be  attributed  to  an  increase  in  x,f.  The  DC  field  range  2  to  3  koe,  which  embraces 
the  spin  flop  value,  was  examined  very  carefully,  but  there  was  no  apparent  instability  effect. 

It  was  estimated  that  the  maximum  intensity  of  the  RF  magnetic  field  in  the  coil  was  0.75  oe, 
which  should  have  been  sufficient  to  exceed  the  instability  threshold  provided  that  was  not 

greater  than  0.375  oe. 

Recognizing  the  probable  difficulty  of  lining  up  the  magnetization  with  the  applied  DC  field, 
as  concluded  from  the  DC  susceptibility  measurements,  the  negative  experimental  result  was 
not  taken  to  imply  that  the  spin  wave  linewidth  was  much  larger  than  expected,  but  that  the  basic 
requirements  for  parallel  pumping  were  not  satisfied. 

The  direction  of  with  respect  to  the  coil  axis  was  changed,  with  the  hope  that  this  would 

cause  M  to  line  up  parallel  to  the  RF  field,  and  the  experiment  was  repeated.  However,  no 
combination  of  magnitude  or  angle  of  applied  field,  or  magnitude  of  RF  field,  resulted  in  the 
observation  of  spin  wave  instabilities. 
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CHAPTER  VI 
CONCLUSIONS 


Magnetic  susceptibility  measurements  on  RbMnF^  samples  yielded  a  value  for  the  exchange 
field  Hpx  in  good  agreement  with,  but  about  7  percent  less  than,  previously  published  data.  The 
measured  value  for  the  Neel  temperature,  94°K,  was  significantly  higher.  The  susceptibility, 
as  a  function  of  applied  DC  field,  showed  no  abrupt  discontinuity,  although  an  analysis  of  the 
static  equilibrium  position  of  the  magnetization  indicates  that  this  is  possible  under  certain 
conditions. 

Antiferromagnetic  resonance  theory,  together  with  the  static  equilibrium  analysis  for  the 
restricted  case  of  lying  in  a  (110)  plane,  has  permitted  prediction  of  the  AFMR  spectrum. 

The  theory  is  in  fair  agreement  with  the  experimental  data  for  the  [111]  and  [100]  directions. 
Observed  resonances  for  the  case  of  Hj^  parallel  to  the  [111]  axis  are  consistent  with  the  in¬ 
terpretation  that  there  were  few,  if  any,  spins  parallel  to  H^.  The  negative  result  of  an  ex¬ 
periment  designed  to  excite  spin  wave  instabilities  by  the  parallel  pumping  technique  is  consistent 
with  this  interpretation. 

The  value  of  the  quantity  JlSH  (H  +  H^),  obtained  by  extrapolating  the  experimental  AFMR 
data  to  zero  DC  field,  is  about  3.5  percent  lower  than  the  value  calculated  from  published  data. 
This  discrepancy  is  consistent  with  that  found  in  the  susceptibility  measurements. 
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APPENDIX 

DERIVATION  OF  THE  TRANSFORMATION  MATRIX  [T] 


Let  (0,X,  Y,  Z)  be  the  coordinate  system  referred  to  the  cube  edges,  and  let  the  (X,  Y,  Z) 
axes  be  rotated  through  an  angle  ©  about  a  line  passing  through  the  origin  and  perpendicular  to 
the  (110)  plane  of  interest,  i.e.,  about  a  [110]  axis  (see  Fig.A-1).  Also,  let  (o,  x,  y,  z)  be  the 

A  A  A  a  A  A 

new  coordinate  system,  and  let  X,  Y,  Z  and  x,  y,  z  be  unit  vectors  along  the  respective  axes. 

A  /\ 

Further,  let  V  be  a  unit  vector  along  the  [110]  direction  in  the  (O,  X,  Y,  Z)  system  with  m  the 

A 

corresponding  unit  vector  in  the  (o,  x,  y,  z)  system,  and  let  U  be  a  unit  vector  perpendicular  to 

A 

the  (110)  plane  at  the  origin  in  the  (O,  X,  Y,  Z)  system  with  t  the  corresponding  unit  vector  in 
the  (o,  x,  y,  z)  system. 

Then 

A  A 

y=  Y  +  X 

V2 


A  ,  A  A  A 

&  =  t  = 

\]~2  VT 

where 

A  A 

U  =  £ 

A  /\  A 

V  =  z  sin  0  +  m  cos  0 

A  a  /\ 

Z  =  z  cos  0  —  m  sin  0 


It  follows  that 
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Y  +  X 
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(X  +  y) 
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Z  =  z  cos  0 


_  (y  +  x) 
\/"2 


cos  0 


sin  0 
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By  rearranging  the  last  three  equations  and  writing  them  in  matrix  form, 
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sin  8 
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Hence, 
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